Codimension-two points in annular electroconvection as a function of aspect ratio 
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We rigorously derive from first principles the generic Landau amplitude equation that describes 
the primary bifurcation in electrically driven convection. Our model accurately represents the 
experimental system: a weakly conducting, submicron thick liquid crystal film suspended between 
concentric circular electrodes and driven by an applied voltage between its inner and outer edges. 
We explicitly calculate the coefficient g of the leading cubic nonlinearity and systematically study 
its dependence on the system's geometrical and material parameters. The radius ratio a quantifies 
the film's geometry while a dimensionless group V, similar to the Prandtl number, fixes the ratio 
of the fluid's electrical and viscous relaxation times. Our calculations show that for fixed a, g is & 
decreasing function of V, as V becomes smaller, and is nearly constant for P > 1. As "P ^ 0, g ^ oo. 
We find that g is a nontrivial and discontinuous function of a. We show that the discontinuities 
occur at codimension-two points that are accessed by varying a. 



I. INTRODUCTION 

With the substantial progress realized in the field of 
pattern formation over the last 15 years has come an 
increased need to make stricter comparisons between 
experiments and their theoretical descriptions Jlj. It 
has become increasingly important for first principles 
theories and numerical models to closely mirror exper- 
imental systems in order to make unambiguous compar- 
isons. A broad range of phenomena is exhibited by lab- 
oratory pattern-forming systems, for example, station- 
ary and traveling patterns, spiral defect chaos, localized 
structures, and so on. An experimental system is typi- 
cally specified by several dimensionless control parame- 
ters which span regions where different patterns are ob- 
served. These regions are bounded by lines or planes in 
parameter space. Codimension-two (CoD2) points oc- 
cur at special values of the control parameters, and are 
the nonequilibrium analogs of multicritical points. Near 
such points, especially interesting and complex pattern 
interactions may be expected 01- 

One of the most successful approaches used to study 
patterns is the Landau amplitude equation for pattern 
amplitudes near a bifurcation. Since the equation can be 
deduced from symmetry, it has found broad applicabil- 
ity in several experimental systems including Rayleigh- 
Benard convection (RBC), Taylor vortex flow (TVF), 
and electrohydrodynamic convection in nematic liquid 
crystals (EHC) 3]. In this paper, we apply the Lan- 
dau amplitude formalism to thin film electroconvection in 
an annular geometry, a system that has previously been 
studied both experimentally^JM 0, S 111 and theoreti- 
cally at the linear stability level^. This system has a rich 
bifurcation behavior, with numerous CoD2 points 0,113. 

The Landau amplitude equation can be rigorously 
derived from the complete set of underlying dynami- 
cal equations by perturbative expansions about the bi- 
furcation point. Very different physical systems with 
a common symmetry-breaking bifurcation have, up to 
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FIG. 1: The coordinate system and film geometry. 



parameter-dependent coefficients, identical amplitude 
equations. Whereas this universal description has been 
tremendously successful in describing the pattern near 
the bifurcation, its quantitative verification relies heavily 
on comparing absolutely the parameter-dependencies of 
the measured and calculated coefficients. This strategy 
has seldom been executed and to our knowledge only in 
RBC, TVF, and EHC has it been generally successful 
3]- Whereas these systems are three-dimensional (3D), 
electroconvection in an annular fluid film, the system dis- 
cussed here, is two-dimensional (2D). RBC, TVF, and 
to a lesser extent EHC, in large part owe their spatio- 
temporal richness to the extra spatial dimension. The 
comparatively simpler spatio-temporal structure of 2D 
thin film electroconvection is accentuated by its special 
geometry. The annular film results in a naturally pe- 
riodic experimental system. This azimuthal periodicity, 
and thus the absence of lateral boundaries, simplifies the 
theoretical treatment and invites interesting experimen- 
tal scenarios such as interposing convection with shear in 
a closed channel d, Q, • 

In this paper, we present a detailed study of the pri- 
mary bifurcation to electroconvection in an annular film. 
The film geometry is shown in Fig. ^ We have previ- 
ously shown experimentally that the system is adequately 
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modeled by a Landau amplitude equation with a cubic 
nonlinearity Q . Our aim here is to rigorously derive the 
parameter-dependence of the coefficient g of the cubic 
nonlinearity and compare it with existing measurements 
from experimental data. In Section II Al we introduce 
the experimental system, its systematically variable di- 
mensionless parameters, and summarize the protocol by 
which g was measured. The theoretical model comprising 
the underlying physics is briefly introduced in SectionllTI 
We proceed to set up the mathematical formalism to re- 
duce the more complicated basic equations to the Landau 
amplitude equation truncated at cubic order in the per- 
turbation expansion. We relegate details of the multiple- 
scales expansion to the Appendix. We study how g varies 
with the two dimensionless parameters and compare the 
trends with experiments in Section IIIII We discuss the 
implications of our work and present a brief conclusion 
m Section im 



A. Previous experimental results 

Electroconvection in an annular film has been the sub- 
ject of several experimental and theoretical studies which 
have examined convection patterns near onset, interac- 
tion with shear flows, linear stability, bifurcations near 
CoD2_points, and more recently turbulent convection 
iillii. 

The electroconvection cell is shown schematically in 
Fig. n It consists of an annulus bordered by two 
concentric stainless steel electrodes with inner (outer) 
radii (ro) ^ 1 cm. A film of smectic A octyl- 
cyanobiphenyl (8CB) liquid crystal doped with tetra- 
cyanoquinodimethane (TCNQ) spans the annulus. The 
resulting film is a weakly conducting 2D annular disk of 
width d = To — Ti, radius ratio a = ri/vo, and thick- 
ness s 0.2 /im. Thickness inhomogeneities relax in the 
freely suspended film which retains its thickness unifor- 
mity even when convecting. The fluid response is New- 
tonian and the material parameters of the liquid crystal 
are well characterized by its 2D mass density p, molecu- 
lar viscosity ?/, and electrical conductivity a. The cell is 
housed in a vacuum chamber which doubles as a Faraday 
cage. An experiment consists of drawing a uniform film, 
placing it under a vacuum, and applying a dc voltage V 
to the inner electrode while holding the outer electrode 
at ground potential. The current / through the film is 
measured. By varying V a current-voltage characteristic 
is obtained. More details regarding the experimental ap- 
paratus and the data acquisition procedure are given in 
Ref. 0. 

A representative I — V curve is shown in Fig. |21 A 
convection threshold at a critical voltage V = Vc \s 
clearly observed. The current is transported by ohmic 
conduction for V < Vc while convection contributes for 
V > Vc- An experimental realization is categorized by 
the dimensionless parameters a, already introduced, and 
'P = eo?y/pcrc?, where eq is the dielectric permitivity of free 
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FIG. 2: (Color online) A representative current-voltage char- 
acteristic (o), the corresponding amplitude-e curve (□), and 
the fit to the Landau amplitude equation (solid line). In this 
case, we have a = 0.47, V = 20.4, and best-fit parameters 
g = 1.9±0.1, h = 5.1±0.2, and / = 0.008±0.002. The dashed 
line is a linear fit to the current- voltage data for V < Vc, in- 
dicating the region of ohmic response. At the onset of elec- 
troconvection, the data depart from the ohmic response. 



space. The parameter V is the ratio of time scales of elec- 
trical and viscous processes of the film and is analogous 
to the familiar Prandtl number of RBC. The experiments 
explored various a by using different combinations of in- 
ner and outer electrodes. At each a, several films of dif- 
ferent s and consequently different P were investigated. 

The I~V data can be expressed in terms of the reduced 
amplitude A and the forcing parameter e. 



1 and e = — 



1 . 



(1) 



where Ir, 



is the contribution to the total current / 



that is due to ohmic conduction. Below the onset of con- 
vection, / = Icond and A = 0. When V > Vc, then 
I > Icond- This results in e > and A > 0. The ampli- 
tude A and control parameter e are then effectively an 
order parameter and a reduced temperature for electro- 
convection. Having transformed the I ~ V data using 
Eq. n the A — e data is then modeled with the phe- 
nomenological steady state Landau amplitude equation 



eA - gA^ - hA^ + f = . 



(2) 



In the modeling, g, h, and / are fit parameters with / 
restricted to be positive. Since g can be negative, the 
quintic term in the amplitude equation is necessary. De- 
tails regarding the data modeling have been reported in 
Ref. 0. 
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Our interest in this paper is to investigate the varia- 
tion of g with a and V. Previous experiments Ifij have 
explored the regimes 0.33 < a < 0.80 and 1 < < 150. 
Measurements of g were found to be roughly independent 
of V except a.t a — 0.33 and 2 < < 8, where g was 
seen to increase with increasing V. As a function of a, 
g was found to be generally increasing with increasing 
a in overall trend. In the next section we recount fea- 
tures of the physical model and set up the mathematical 
formalism for perturbatively solving the equations. 



II. MATHEMATICAL FORMALISM 

A linear instability mechanism gives rise to electro- 
convection in a freely suspended fluid film 0, 0|. The 
film, an annular disk on the z = plane, is subjected 
to electric potential boundary conditions of V volts at 
its inner edge (r = r^) and zero volts at its outer edge 
(r — To). The inner electrode r < is at F volts and the 
outer electrode r > Tq is at zero volts. The potential is 
zero at infinity. The electrostatic boundary value prob- 
lem prescribed by these conditions, with the film being a 
conducting liquid, implies that a surface charge density 
develops on the film's upper and lower free surfaces. Pos- 
itive charge preferentially accumulates close to the inner 
positive electrode and negative charge at the grounded 
outer electrode. Consequently, an electric force is ex- 
erted on the fluid by the action of the radially outward 
component of the electric field on the surface charge den- 
sity. Note the striking analogy with RBC where thermal 
boundary conditions on a thermally conducting liquid 
lead to hotter, less dense fluid near the hot boundary and 
colder, more dense fluid near the cold boundary. The ac- 
tion of the gravitational field on the mass density exerts 
a force on the fluid. In both cases when the forcing over- 
comes dissipation, the fluid becomes linearly unstable to 
perturbations, resulting in convection. 

The relationship between the electric potential and the 
surface charge density at any spatial position on the film 
is nonlocal in electroconvection, unlike thermal convec- 
tion where the temperature and mass density are locally 
related. The surface charge density is directly related to 
the discontinuity in the component of the electric field 
perpendicular to the film at z = 0. As a result, electric 
fields in the full 3D space determine the charge density 
of the film at z = 0. 

We use the cylindrical coordinate system (r, 9, z). The 
film is in the z = plane with radial coordinates between 
Ti < r < To and has the areal material parameters de- 
scribed in Sec. II Al It is assumed that the fiuid is uniform 
in temperature and there is no significant ohmic heating. 
Since the velocity field is 2D as the film fiows in the z = 
plane, we choose to describe it implicitly with the stream- 
function ipir^O). The surface charge density is q{r,9) 
and the electric potential is denoted by ^'3(r, 6*, z). The 
electric potential in the film is for convenience denoted 
^{r,9) = ^3{r,9,z = 0). The fluid is described by the 



incompressible Navier-Stokes equation with an electrical 
body force. The surface charge is adverted by the flow 
and transported by ohmic conduction. Laplace's equa- 
tion is obeyed by the electric potential for z 7^ 0. Nondi- 
mensionalizing lengths by the fllm width d ~ ro ~ ri, 
electric potential by V, and time by eod/a in the mo- 
mentum, charge conservation, and Laplace equations we 
get the set of governing equations 





Vdt 

1 



(V X V X 0) + 7^ ^V* X Vq^ 



^^(Vx0).vj(Vx Vx 



dq ,^ 



Vg- V''^' = 0, 



V^^-s = 0, 



q = -2- 



dz 



(3) 

(4) 

(5) 
(6) 



z=0+ 



where the dimensionless parameters are 



7^ ; 



ar/ 



and 



£0^7 
pad 



(7) 



Details regarding the assumptions and the derivation can 
be found in Refs. HQ. 

In the standard manner 0], we write the streamfunc- 
tion, charge density, and electric potential as the sum of 
the base state solution and a perturbation. See Ref. 
for details. The resulting equations can be succinctly 
expressed as 



£.C = B 



Here we have the operator 



C = 



( V4 



V 




1 









(8) 






^3 



(9) 

We have denoted functions of the base state by super- 
script (0). Note that = drr + + ^dee while 
V3 = + dzz- The functions C and B are 



C = 



0(r,0) 
q{r,9) 
*(r, 9) 
*3(r,0,z) 



(10) 



and 

^= dtq+l[{drq){dgcl,)~{dgq){drCl))] 



V / 

(11) 

The multiple-scales perturbation theory employed in 
our treatment is the same as that given in Ref. 11] for 
electroconvection in a rectangular geometry. We expand 
Eq. fusing the slow time scale T = et, where e — TZ/TZc — 
1 is the reduced control parameter defined earlier in Eq.^ 
Collecting terms of the same expansion order in e we 
write 

^qCq = Bq , 
CqCi + CiCq = Bi , 
C0C2 + CiCi + C2C0 = B2, (12) 

for orders e^/^, e, and e"^/^. A systematic sequential solu- 
tion of the above equations results in the necessary con- 
dition 

FidTA + F2A + FiAlAl"^ ^0. (13) 

The amplitude equation in the fast variables is of the 
Landau form 

TdtA ^ eA- gA\A\^ , (14) 

where 

r = and g = . (15) 

i^2 ^2 

The functions Fi, F2, and F4 are further discussed in the 
Appendix. Other details regarding the numerical evalu- 
ation of g and r are also given in the Appendix. 

III. DISCUSSION 

In this section we present the results of our calculations 
of the functional dependencies of g = g(a,V). Since the 
experiments on electroconvection are, due to the values of 
the physical and dimensional parameters, constrained to 
large V > \, we start by discussing g = g{a) for large and 
essentially infinite V. In figureOlis plotted the calculated 
values of g for 0.6 < a < 0.8. In overall trend, the Lan- 
dau cubic coefficient g increases with a and approaches a 
limiting value as a — > 1. At a = 0.8, g = 2.570 is within 
10% of the limiting value g — 2.842 calculated earlier 
in a Cartesian or rectangular geometry. The coefficient 
g(a) decreases with a for intervals over which the criti- 
cal mode m = is constant. This trend is punctuated 




a 

FIG. 3: g versus a for P = 123. 

by discontinuities in g at values of a where both m and 
m4-l are equally unstable; these are the codimension-two 
(CoD2) points. At these jumps, the value of g increases 
and more than compensates for the region of decrease at 
each m. The result is an overall increasing trend. 

TableQlcompares experimental measurements of g from 
Ref. 0] at six radius ratios with the results of the present 
calculations. First note that in spite of the large scatter 
in the experimental measurements, they still show the 
overall increasing trend of g with a. Further, the exper- 
iments have widely separated values of a. and so are not 
able to resolve the discontinuities in g for which further 
experimental work in a restricted but densely sampled 
range of a would be needed. And finally, in comparing 
measurements to the calculations, there is a significant 
disparity on the order of 10 — 30%. Interestingly, the ex- 
perimental measurements show that as the radius ratio 
is decreased, the Landau cubic coefficient becomes neg- 
ative. This implies that there is at least one tricritical 
point [g — 0) that demarcates the super- and subcriti- 
cal branches. The current calculations are cumbersome 
to extend to smaller a because this requires additional 
orders in the expansion of the stream function |12|. 

The CoD2 discontinuities in g are larger at smaller 
a as is evident from Fig. |31 To further quantify this 
observation, we have graphed the fractional discontinuity 
in g, 2(5,„+i - gm)/{gm+i + 5m), versus the fractional 
change in the mode number, 2(m -t- 1 — m)/{m -I- 1 -I- 
m), at several CoD2 points, as shown in Fig. 01 In the 
above, gm is the value of g for a particular critical mode 
TO — TO^. For the CoD2 point at m = 8,9, the jump 
in g is about 3% suggesting that experiments looking 
at the behavior near CoD2 points will have typically to 
measure g with 1% resolution. The current resolution of 
experimental measurements is about 10% and will have 
to be significantly improved if any meaningful study is 
to be made. On the other hand, since the fractional 



a 


g (experiment) 


P 


g (theory) 


0.33 


-0.74 ± 0.23 


2.1 < P < 4.4 




0.47 


1.64 ± 0.06 


13.5 <r < 20.7 




0.56 


0.73 ± 0.15 


59.4 <V < 100.8 




0.60 


2.72 ± 0.34 


31.3 <r < 38.9 


2.372 


0.64 


1.87 ± 0.10 


25.2 <r < 63.0 


2.417 


0.80 


2.21 ± 0.29 


15.3 <V < 142.8 


2.570 


1.00 ('plate') 




P = oo 


2.842 



TABLE I: Experimental measurements of the coefficient of 
the cubic nonlinearity, g. The theoretical g at a = 0.60, 0.64, 
and 0.80 are for V = 123. 



discontinuity in g increases strongly with the fractional 
discontinuity in mode number, it is conceivable that at 
small mode CoD2 points, say to = 3, 4, the discontinuity 
in g will be very large and measurements with 5 — 10% 
accuracy may suffice. This will involve working at small 
a which is a challenge experimentally since broad films 
will have to be drawn. 

We now discuss the T'-dependence of g with constant 
a. The calculated values of 5 for 0.001 < V < 10 for 
a = 0.674 and 0.676 are plotted in Fig. 13 These a 
straddle the CoD2 point for to = 10, 11. The cubic Lan- 
dau coefficient is practically constant for V > 0.1. How- 
ever, for small V, g diverges as P — > 0. The inset in 
Fig. 13 shows the absolute value of dg / dV as a function 
of V. As V decreases, the discontinuity in g at the CoD2 
point TO = 10, 11 increases as is seen by the two diverging 
curves in Fig. |21 For V > 1, the fractional discontinu- 
ity is about 2% growing to 10% for V = 0.01 and 30% 
for V — 0.001. This suggests that experimental measure- 
ments near CoD2 points are better performed at small V. 
The combination of small a and V are required for large 
changes in g through CoD2 points. Unfortunately, exper- 
iments can seldom access this regime of parameter space. 
Most experiments have been performed for a > 0.30 and 
V > 1. The latter restriction is primarily because the 
material parameters are only slightly adjustable. 

Langford and Rusu have considered several qualitative 
bifurcation scenarios for this system They focused 
on the possible transitions at CoD2 points. Our work 
indicates that the bifurcations at CoD2 points are su- 
percritical for most radius ratios. Subcritical transitions 
may be possible at small a [T^ . 



IV. CONCLUSION 

In summary, we have examined the variation of g, the 
coefficient of the nonlinear term in the amplitude equa- 
tion Eq. 1141 with radius ratio a and Prandtl number V 
for electroconvection in a two-dimensional annular fluid. 
We have observed that the steady-state amplitude of 
the bifurcation from a conducting to a convecting fluid 
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FIG. 5: (Color online) The coefficient g versus P at a = 0.674 
(solid line) and a = 0.676 (dashed line). The inset shows 
\dg/d'P\ versus V for the same a. 



can change discontinuously as a passes through CoD2 
points. The dependence of g on P is not strong, except 
for V < 0.1. 

In the experiments considered here, the variation of g 
is a consequence of changing the geometry of the system 
via the radius ratio a. It is also possible to pass through 
CoD2 points by applying a Couette shear to the fluid 
at fixed a J]. Experiments revealed that g was found 
generally to decrease with increasing shear and discon- 
tinuities in g were observed at CoD2 points Apply- 
ing the methods employed in this paper to the case of 
nonzero shear, it should be possible to establish whether 
codimension-three (CoD3) points exist at special values 
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of the parameters [ij 
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APPENDIX: MULTIPLE-SCALES EXPANSION 

We used a multiple-scales formalism analogous to that 
developed for RBC in Ref. 0, in which a compatabil- 
ity condition is derived for the amplitude of the slowly 
varying envelope of the convection rolls. We generally 
adhere to the notational conventions of Ref. and our 
treatment is similar to that of Ref. ) adapted to the 
annular geometry and including the effects of finite V. 
The multiple-scales expansion of Eq.|Slgives Ea.ll4lwhere 
e - (7^ - 7^e)/7^c, r = -F1/F2, g = ~Fi/F2, 



F, ^ ^m°7^°(^[5.vI/("'go-a,.<z(°)'^o]) 



(A.l) 
(A.2) 



im 



0/gfco"^c I 



2dr qQ(t)\ + q^drCpl 
-2qldAl]) , 

.) = r rdr (...). 



(A.3) 
(A.4) 



Complex conjugation is represented by a superscript *. 
The functions that appear in Eas. lA.llA.3l are as follows. 
The solutions of the linear stability problem at to = toJ^ 
are 



aim'ie 



^0 = Ep ^p0op 90 = Ep ^p9op ^0 = Ep -^p^-op 



;p 7 



lop = iy^vqpiqmo-i 
I 

*0p = »y^^OpiV'm»;l 



(A.5) 
(A.6) 
(A.7) 

(A.8) 



The expansion functions Cm;p(r) and il>m;i{T,z) satisfy 
the boundary conditions in polar coordinates Q and 

gm;/(7') = -2[a^Vm;/(7', z)]|2=o+- Thc fuUCtioUS 

01 = 92 = E; bilo-i *2 = E; bitpO;i 

(A.9) 



p 
,m 



'^mpl 



drq^°^Cm:.p'>Pra-l) , 



N2 = 
N3 = 



mini 



(qadrlia - drqa^o) 



rV 







^-^[drqocpo - qodrM , 
r 

ITJ1P — — 

-[{drqo<l>0 - drqo(t)o) 

r 

~{qodr4>Q - qo9r4>o)] , 



(A.ll) 
(A.12) 

(A.13) 
(A. 14) 

(A.15) 



satisfy the order e multiple-scales equations. The coeffi- 
cients Ep are specified by 

'^TkpEp = {NiC2mO;k) 

p 

-2TO°7^° ^(iV2V'2m«;;)2'2m0fc;/X2m0/ , 

(A.16) 



T, 



kp 



P2m0p^kp 



S2m°pl Z2ni°kl/Xlm°l >(A.17) 



; Cm: 



Zmkl - (^(a,.Vl'(")q™;Z-9,q(°Vm;/)) .(A.18) 
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The solutions of the adjoint equations evaluated at m 



where 



ho — (pboe <= Qbo — Qboe <= ^to — ^boe •= 
ho = J2p Bpcjjbop qbo — Yl,p BpQbop ^'ho ~ J2p Bp^bop 

(A.19) 

qbOp = ^mO;p , (A. 20) 

(j^bOp = iy^^VbOplCmO-i , (A. 21) 



VbOpl = -Sm°lp/Ptn«l , 



"ifbOp = 



The coefheients Bp are derived from 



The amplitude A is normalized by setting 

{{qur}) 



Nu- 1 



((rf))) 



(A.22) 



-m°n%¥°^i,Op ■ (A.23) 



J2Bp[--m''^n%q^°^$bop+^^qbop+2{d,^Mp)\.=o+] = 

(A.24) 



(A.25) 



Ztt 



2tt 



dO I rdr{. 
h) 



(A.26) 



is the radial component of the velocity field, 

a is the conductivity, and is the radial component 
of the base state electric field. The Nusselt number Nu 
is the ratio of the total current density to the conducted 
current density, spatially averaged 



sity 

[nj. 



To find the coefficients of the normalized amplitude 
equation, Eg. 1141 we evaluate Eas. lA.llA.3l for a given a 
and V. The charge density expansion functions g™;/ are 
computed via an approximation of the electrostatic equa- 
tions in which 5" and q are linearly related ■ The series 
solutions of the linear (Eq. IA.5|) and adjoint (Eq. IA.19|) 
problems are terminated at p = 1. The number of ex- 
pansion functions employed in the linear fEas. IA.6IA.8p . 
first-order (Eq.EH), and adjoint fEas. IA.20lO3)l solu- 
tions is twenty. 



